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SrR.U 2 06 has unique magnetic properties. It is characterized by a very high Neel temperature, 
despite its quasi-two-dimensional structure, and has a magnetic moment more than twice reduced 
compared to the formal ionic count. First principles calculations show that only an ideal Neel 
ordering in the Ru plane is possible, with no other metastable magnetic solutions, and, highly 
unusually, yield dielectric gaps for both antiferromagnetic and nonmagnetic states. We demonstrate 
that this strange behavior is the result of the formation of very specific electronic objects, recently 
suggested for a geometrically similar Na 2 lr 03 compound, whereby each electron is well localized 
on a particular Ru6 hexagon, and completely delocalized over the corresponding six Ru sites, thus 
making the compound both strongly localized and highly itinerant. 


The recently discovered PQ SrRu 2 0 6 has attracted con¬ 
siderable attention because, despite being a very 2D ma¬ 
terial, it shows an exceptionally high Neel temperature 
of -560 K |2 E]. As we will argue in this paper, this 
is by far not the only, and maybe not even the most 
intriguing property of this material. Ru 5+ has a half- 
filled t 2 g electronic shell, and exhibits insulating behav¬ 
ior. Naturally, it was interpreted as a Slater insulator 
(maybe Mott-enhanced), with Ru in the high spin state, 
S = 3/2. However, the experimentally measured ordered 
magnetic moment is only M = 1.3 — 1.4 ^b ED®, 2.3 
times smaller than expected for S = 3/2 (M = 3 fi b )• 
This was ascribed to hybridization with oxygen [2H1], but 
it should be noted that such strong suppression of mag¬ 
netic moment in a good insulator is unheard of. Even in 
the metallic SrRuCL the hybridization suppresses the to¬ 
tal magnetic moment of Ru 4+ only from 2 to 1.7 [Ib, and 
in Sr 2 YRuC>6 Ru 5+ has essentially exactly 3 /is, with ba¬ 
sically the same Ru-0 distances as in SrRu 2 C>6 [Ej- Hiley 
et al [2| mention the case of LLRuCL l6j, where a sup¬ 
pression down to M = 2.0 hb was reported for the same 
oxidation state, which is, however, still twice smaller a 
reduction compared to SrRu 2 06, and the material might 
actually be a metal (no transport data have been pub¬ 
lished) . 

Electronic structure calculations HEj so far have not 
resolved the mystery, but have only added to the con¬ 
fusion. It was found that only the ideal Neel state can 
be stabilized in the calculations, even though ions with 
S = 3/2 are usually very stable and while they disor¬ 
der with temperature, never lose their magnetic moment 
completely. At the same time the moment found in the 
calculations matches the experimentally measured one 
within 8%, suggesting that the role of Coulomb corre¬ 
lations beyond the standard density functional theory 
(DFT) is negligible [7]. The instability of the ferromag¬ 
netic (FM) state was traced down to the presence of a 
dielectric gap in nonmagnetic calculations [I], but that es¬ 


sentially translates one mystery into another: why does 
a highly symmetric Ru sublattice, with no dimerization 
or clusterization, with a half-filled t 2 g band, show a size¬ 
able nonmagnetic gap? Singh mentions [4] that the gap 
is allowed by symmetry, since the unit cell includes two 
Ru atoms that can, in principle, form a bonding and an 
antibonding band, but does not elaborate about how a 
structure with each Ru having three equivalent bonds 
manages to develop a bonding-antibonding splitting. 

Similarly, it was pointed out that, even though 
SrRu 2 C>6 is extremely 2D magnetically, there is still some 
residual interlayer coupling, J±AI 2 ss 1.5 meV, as well as 
a single-ion magnetic anisotropy, estimated to be ss 1.4 
meV/Ru [2]. It was suggested that the anisotropy ]3| or 
interlayer coupling [3] are responsible for the large Tjv, 
implying that the (unknown) mean field transition tem¬ 
perature is extremely high. Tian et aZ[3] attempted to 
describe this system by a three nearest neighbor Heisen¬ 
berg model with parameters derived within the pertur¬ 
bation theory in the limit of a Hubbard U much larger 
than the hopping, U t. However, the fact that fer¬ 
romagnetic arrangement is completely unstable (in fact, 
as we show below, no parallel nearest neighbor moments 
are stable), indicates that the system is strongly non- 
Heisenberg, casting very strong doubt on the relevance 
of such models. Additionally, the fact that the system 
is very weakly correlated makes such a perturbation the¬ 
ory unphysical. Hiley et al [2] used a hybrid functional 
that overestimates the equilibrium magnetic moment and 
thus the exchange parameters, [5] as well as yields a very 
large band gap of 2.15 eV, totally inconsistent with the 
observed weak temperature dependence of the resistivity. 

An explanation of all these oddities can be consistently 
found in the so-called molecular orbitals (MO) picture, 
which was first brought up in connection to Na 2 IrO 3 0 
and later found also in RuCL pH]. Basically, this picture 
is based on the idea that for ideal 90° Ru-O-Ru bond 
angles (the actual angles are 101°) the O-assisted Ru-Ru 
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FIG. 1. Density of states (DOS) projected on molecular 
orbitals of different symmetries in nonmagnetic GGA calcu¬ 
lations (WIEN2k results). The Fermi energy is set to zero. 

hopping is only allowed for one particular pair of the t 2 g 
orbitals for each hexagonal bond, denoted t\ in Ref. [5]. 
If all other hoppings are neglected, it leads to a curious 
situation where every electronic state is fully delocalized 
over a particular hexagon, but never leaves this hexagon. 
One can say that the electrons are fully localized (form 
nondispersive levels) and fully delocalized (each state is 
an equal weight combination of six orbitals belonging to 
six different sites). If a direct overlap of the t 2 g orbitals 
(which always exists in the common edge geometry) is 
included, as well as deviations of the angle from 90°, two 
more hoppings emerge: one between the same orbitals 
on the neighboring sites, ti, and the other an O-assisted 
second neighbor hopping between unlike orbitals, t 2 . As 
long as t\ is dominant, the MO model still applies, and 
can be readily solved. The solution entails six bands, A lg , 
E 2u , Ei g and B iu (the E bands being double degenerate 
in each spin channel), whose dispersion is controlled by 
ti, and whose centers are located at 2 (t[ + t' 2 ), (t[ — 
t 2 ), —(ti + t 2 ) and —2 (t[ — t 2 ), respectively. In Na 2 lr 03 
t\ ss — 3t 2 , so that the A\ g and E 2u practically merge. 
This accidental degeneracy also leads to much stronger 
spin-orbit effects than would have been possible had the 
MO bands remained well separated, and to considerable 
destruction of the MO picture in the relativistic case. On 
the other hand, the hopping parameters for SrRuoOe, as 
calculated in Ref. m. are similar to those in Na 2 lr 03 , 
in the sense that again t\ = 300 meV is by far the largest 
hopping, and the only other sizeable hoppings are t\ = 
160 meV and — 1 2 ~ 100 — 110 meV. Note that here \t 2 \ 
is again about 1/3 of t\. Thus, the A lg and E 2u bands 
merge, while E- tg and Bi u remain separated, as one can 
see in Fig. [T] Projecting the density of states onto MOs, 
we observe that the predicted characters are very well 
reproduced. The distance between the centers of the E 2u 
and E\ g bands is about 0.8 eV, and their width is about 
0.6-0.7 eV, thus providing for a small gap of « 50 meV. 

It is instructive to compare SrRu 2 Og with Na 2 lr 03 
and with L^RuOs. All these compounds share the same 
crystallographic motif, but feature a different number of 


d electrons: 5, 4 or 3. In the iridate, a single hole in the 
upper Aig singlet is prone to both strong correlations 
and, due to near degeneracy between Ai g and E 2u , to 
spin-orbit interaction. As a result, as one increases the 
spin-orbit coupling, the A\ g singlet is gradually trans¬ 
formed into the j e f / = 1/2 singlet [9]. Either way, a 
half-filled singlet triggers Mott physics even if the Hub¬ 
bard U is small. This transformation controls most of 
the interesting physics in this compound. Li 2 Ru 03 has 
two d holes, providing it with an opportunity to form 
strongly bound covalent dimers. This is exactly what 
happens, and the MO on the hexagons transforms to an 
MO on the Ru dimers resulting in the spin singlet ground 
state .IS]. Neither Mott nor spin-orbit physics is rele¬ 
vant on the background of the strong covalent bonding 
in dimers. Finally, SrRi^Og has the six MO bands half- 
filled, and the gap is formed between the lower and the 
upper MO triads. Similar to L^RuOs, both Mott and 
spin-orbit effects are of minor importance, and the gap 
structure inherent to the MO picture gives rise to unique 
magnetic properties. 

Let us now turn to the energetics of the material. First, 
we have confirmed, using the WIEN2k package nans, 
the numbers published by SinghjTj regarding the inter- 
planar coupling, single-site anisotropy and Ru magnetic 
moment. We also confirmed that the ferromagnetic struc¬ 
ture cannot be stabilized. Moreover, the so-called stripy 
and zigzag magnetic patterns [9], where one or two out 
of three bonds are ferromagnetic, and the net moment 
is zero, cannot be stabilized. This indicates that besides 
the obvious influence of the nonmagnetic gap there are 
other factors strongly disfavoring ferromagnetic bonds. 
In fact, given that the gap value is ten times smaller than 
Ru Stoner factor [5], and the calculated magnetic moment 
in the Neel state is ~1.3 it is surprising that the fer¬ 
romagnetic bonds do not stabilize with a finite moment. 

In order to gain more insight into the problem, we 
turned to the VASP code m ns, which is faster and 
has the capability to restrict magnetic moments to a cer¬ 
tain direction, or to both a direction and a magnitude (we 
confirmed that the energies of collinear magnetic states 
agree with those found in WIEN2k). First, we computed 
the total energy for a canted antiferromagnet (AFM), re¬ 
stricting the angle with the z axis to be ±0 for the two 
Ru’s in the cell. The results are shown in Fig. [2] Note 
that for the largest canting angle we were able to con¬ 
verge, 35°, the energy of the magnetic state is already 
higher than that of the nonmagnetic one. Also note how 
soft the magnetic moments are: despite the sizeable equi¬ 
librium moment, the energy cost of total suppression of 
magnetism is less than 80 meV, only 50% larger than the 
transition temperature. This is, again, an indication of 
the great role of itinerancy, and specifically, delocaliza¬ 
tion over Rug hexagons. 

Interestingly, suppression of magnetism with canting 
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FIG. 2. Magnetic energy (squares) and magnetic moments 
(circles) as a function of the canting angle of spins, starting 
from the Neel antiferromagnetic structure. Results are from 
VASP calculations. 


cannot be described by a naive combination of a local 
Hamiltonian for itinerant magnets[l6j, E = ^2 OiM 2 ®, 

i> 0 

where M is the magnetization, with a Heisenberg term. 
While the total energies at a fixed canting angle <f> < 35° 
can be very well described by this Hamiltonian with just 
three terms, E(M) — E( 0) = a\M 2 + CI 2 M 4 + a^M 6 , not 
only does the first coefficient show appreciable angular 
dependence (as in the Heisenberg model), but also the 
second, and, to a lesser degree, the third. Instead, a 
good fit could be obtained with the following formula: 

E = —81.3 M 2 + 16.9M 4 + 2.0M 6 + 359.2 M 2 sin 2 (</>) 

— 165.8M 4 sin 2 (</>) + 27.6 M 6 sin 2 (</>), (1) 

in meV/Ru. Note that the angle between the moments is 
8 = 7T —2 tj>, and that there are 1.5 times more bonds than 
sites. Thus, the proposed Hamiltonian looks as follows: 

H = {98.3 M 2 - 66.0M 4 + 15.8M 6 } 

sites 

+ {179.6(M • M') - 82.9|M||M'|(M • M') (2) 

n.n. 

bonds 

+ 13.8|M| 2 |M'| 2 (M • M')}. 

The Heisenberg term is extremely strong (JM 2 = 
dH/dcos (9) ss 1600 K), and, without it, local magnetic 
moments fail to form. 

To this Hamiltonian one needs to add a small interlayer 
term ^ , where i and i' belong to the neighbor¬ 

ing planes, and the magnetic anisotropy J2 -DM 2 , where 
J 1 ss 0.9 meV, and D ss 0.8 meV. 

In principle, at this point one would need to perform 
a Monte Carlo simulation using this Hamiltonian and 
determine the transition temperature. However, it is 


notoriously difficult to distinguish a Kosterlitz-Thouless 
phase in a quasi-2D system from the true long range or¬ 
der, so that one should be very skeptical of any Monte 
Carlo simulation that claims to establish a Neel temper¬ 
ature, T/v, without first showing that in the isotropic 2D 
limit Tn truly vanishes. The softness of the moment, 
expressed via Eq. [2j additionally complicates the simu¬ 
lation. We leave this daunting task to more experienced 
Monte Carlo simulators, but mention that the numbers 
that we have deduced are in the right ballpark. For in¬ 
stance, Costa and Pires showed [T? that for the square 
lattice Tn/Tmf ~ 0.8 (D/J) 0 2 . For three neighbors, the 
mean field transition temperature Tmf ~ JM 2 ss 1600 
K, which together with DM 2 « 1.4 meV results in 
Tjv ~ 500 K. On the other hand, for the cubic quasi-2D 
model with J± ^ 0, D = 0, Yasuda et al ng found that 
Tat « 4.27JM 2 /[3.12+log(J/Jj_)], which for our parame¬ 
ters translates into 900 K. Thus, we conclude that (a) the 
Mermin-Wagner theorem is mainly lifted via the interpla- 
nar coupling[3j, and not via the single site anisotropy [1], 
and (b) the softness of the magnetic moment, i.e., longi¬ 
tudinal fluctuations, plays an important role, suppressing 
T/v by up to a factor of two. 

Let us now discuss how and why MOs support a Neel 
antiferromagnetism in SrRu 2 06- In the nonmagnetic 
state, the three lower MO bands, Bi u and Ei g , are fully 
occupied. Imposing uniform spin polarization does not 
change the occupancy of these states, unless the induced 
exchange splitting is larger than the gap, and this is why 
the ferromagnetic order is unstable. On the contrary, im¬ 
posing the staggered magnetic field of ±A does not break 
the MO band structure, but rather increases the gap be¬ 
tween E\ g and E 2 U (hr the lowest order in A, by A 2 /t' 1 ). 
In the same order, we can calculate the change of the oc¬ 
cupancies and find that the spin-up sites acquire magne¬ 
tization of 5A/2t , 1 fiBi and the spin-down sites — 5A/2t , 1 
[iB- The signs are consistent with the assumed signs of 
A, which tells us that with sufficiently large Hund’s rule 
coupling the system will become unstable against such 
a staggered magnetization (but will resist any ferromag¬ 
netic component); of course, quantitative analysis is im¬ 
possible on this level of simplification. Obviously, the 
equilibrium moment can be anything between 0 and 3 
li b ■ It is not “suppressed” from the putative S = 3/2 
state, but is set by the interplay between the Hund’s rule 
coupling on Ru and the details of the density of states 
of MOs. A corollary from the above arguments is that 
the dielectric gap depends quadratically on the Ru mo¬ 
ment; Fig. [3] illustrates that this is indeed the case, to a 
reasonable accuracy. 

Let us emphasize, that the molecular orbitals are not 
just another way to describe the electronic structure of 
SrR.U 2 06, but have profound physical meaning. It is in¬ 
structive to compare it with another recently investigated 
high-T/v material, SrTcOa, where the transition metal 
also has a 4 d 3 configuration and S = 3/2. It was argued 
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FIG. 3. Total density of states (DOS) calculated for several 
values of Ru moments, M, in the fixed-spin-moment proce¬ 
dure for the Neel AFM. The black lines illustrate the fact that 
he band gap is approximately quadratic in M. Results of the 
LAPW calculations. 

m that Tjv is so high because SrTcCU is in an intermedi¬ 
ate regime between itinerancy and localization, which is 
optimal for magnetic interactions. Indeed, LDA+DMFT 
calculations, well suited to this regime, have been per¬ 
formed by Mravlje et al. jT5j, who found T N ss 2200 
K. The experimental number is about 1100 K. To com¬ 
pare this result with SrRu 2 C> 6 , we have also performed 
LDA+DMFT calculations with the AMULET code SO], 
using an effective Hamiltonian constructed for Ru t- 2 g or¬ 
bitals and interaction parameters U = 2.7 and J = 0.3 
eV as calculated in Ref. [3] (parameters for Tc are very 
similar). The corresponding temperature dependence of 
the magnetic moment is shown in Fig. [4j Not surpris¬ 
ingly, we found about the same Neel temperature (2000 
K) as Mravlje et al. [T5], and an even larger magnetic 
moment (M «2.7 vs. 2.5 /is). The difference, however, 
is that experimentally in SrRu 2 C >6 both T)v and M are 
about twice smaller than in SrTcOa. Mravlje et al. as¬ 
cribed their overestimation of TN to nonlocal fluctua¬ 
tions, missing in the DMFT, but observed no reduction 
in the ordered moment at all, while in our case the re¬ 
duction in both T]y and M 2 is of the same order, about 
a factor of 4. This clearly indicates that there is a fun¬ 
damental difference between the two compounds, going 
much beyond just the difference in dimensionality, which 
is related to the presence of MOs in one and their ab¬ 
sence in the other compound. A proper account of the 
molecular orbitals within DMFT can only be done in the 
cluster extension of this method [21], which could shed 
more light on this compound. 

Another interesting question that arises in connection 
with this material is what would happen if it were doped 
with, for instance, a rare earth element. To address this 
scenario, we simulated doping by adding electrons to the 



FIG. 4. Magnetic moment as calculated in the LDA+DMFT 
approach. The continues-time quantum Monte-Carlo (CT- 
QMC) solver [22] was used in these calculations. 

system (with a compensating constant background). The 
energy difference between the FM and Neel AFM states 
decreases upon electron doping, as seen from Fig. [5] 
The FM configuration immediately becomes metastable, 
whereby all doped electrons go into one spin subband, 
rendering the material half-metallic. The ground state 
remains antiferromagnetic, but its energy advantage is 
gradually decreasing. Thus, one expects that the crit¬ 
ical angle <f> (which was ~35° in undopped case) will 
grow with doping, and the Hamiltonian Q will be corre¬ 
spondingly modified; this may result in a rapid change of 
magnetic properties with doping, which deserves further 
theoretical and experimental investigation. 

To summarize, we have found that: 

(i) The electronic structure of SrRi^Og is dominated 
by molecular orbitals. Each electron is, to a good approx¬ 
imation, localized on a particular Rug hexagon, and com¬ 
pletely delocalized over the corresponding six Ru sites. 

(ii) This structure sports an excitation gap that pre¬ 
vents formation of ferromagnetic bonds, but is consistent 
with nearest neighbor antiferromagnetism. The corre¬ 
sponding magnetic interactions cannot be mapped onto 
a localized spin model, be it Heisenberg or biquadratic 
Hamiltonian with arbitrary long range. Neither can it be 
described as purely itinerant magnetism, but features in¬ 
teresting elements of both. This duality reflects the dual 
character of the electronic structure, where electrons are 
simultaneously completely delocalized and strongly local¬ 
ized on the Ru hexagons. A corollary is that any devia¬ 
tion from the collinear Neel order is severely punished by 
kinetic energy, which, in turn, provides for the anoma¬ 
lously large transition temperature. 

(iii) The gaps in the non-magnetic and antiferromag¬ 
netic states have the same nature, and one is continu- 
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FIG. 5. Electron doping dependence of magnetic moments 
(per Ru and total) and total energy difference between FM 
and Neel AFM states on the electron doping. 

ously transformed into the other as the magnetization 
increases. On the contrary, the ionic picture assigning 
the moment of 3 hb to each Ru and associating the 
gap in magnetic states with spin-up/spin-down splitting 
is qualitatively incorrect. The observed and calculated 
magnetic moment of 1.3 /r b is a manifestation of the 
molecular orbital nature of electronic states, and should 
not be viewed as a spin S = 3/2 reduced by hybridiza¬ 
tion. 

(iii) The magnetic properties of doped SrR^Og (e.g. 
by Na or La) are expected to be very different from the 
stoichiometric case. One may anticipate interesting and 
very different physics emerging, which can be a subject 
of forthcoming research. 
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